Spin Liquid Behavior in Electronic Griffiths Phases 
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We examine the interplay of the Kondo effect and the RKKY interactions in electronic Griffiths 
phases using extended dynamical mean-field theory methods. We find that sub-Ohmic dissipation is 
generated for sufficiently strong disorder, leading to the suppression of Kondo screening on a finite 
fraction of spins, and giving rise to universal spin-liquid behavior. 
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Disorder induced non- Fermi liquid (NFL) behavior has 
remained an important focus of study in heavy fermion 
systems According to the "magnetic Griffiths phase" 
scenario |2|, this behavior reflects the formation of rare 
magnetic clusters with large susceptibilities, similarly as 
in disordered insulating magnets |3]- Another approach 
focused on the interplay of disorder and the Kondo ef- 
fect far away from any magnetic ordering 4]. More re- 
cent work 0, @ has demonstrated that such "electronic 
Griffiths phases" are a generic feature of strongly corre- 
lated electronic systems with disorder. Neither picture, 
however, seems satisfactory for the following key reason: 
in both scenarios, the resulting NFL behavior is char- 
acterized by power law anomalies, with non-universal, 
rapidly varying powers. In contrast, most experimental 
data seem to show reasonably weak anomalies, close to 
marginal Fermi liquid behavior [lj. 

Physically, it is clear what is missing from the the- 
ory. Similarly as magnetic Griffithsphases, the electronic 
Griffiths phase is characterized |j,Lj,|l| by a broad distri- 
bution P{T K ) ~ (TV) of local energy scales (Kondo 
temperatures), with the exponent a ~ W~ 2 rapidly de- 
creasing with disorder W. At any given temperature, 
the local moments with Tx{i) < T remain unscreened. 
As disorder increases, the number of such unscreened 
spins rapidly proliferates. Within the existing theory 
0- 0- El these unscreened spins act essentially as free 
local moments and provide a very large contribution to 
the thermodynamic response. In a more realistic descrip- 
tion, however, even the Kondo-unscreened spins are not 
completely free, since the metallic host generates long- 
ranged Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
actions even between relatively distant spins. As the 
RKKY interaction has an oscillatory character, the ef- 
fective interactions between the (randomly located) un- 
screened spins will be random in magnitude and sign. In 
this Letter, we use an extended dynamical mean field for- 
mulation to examine the role of such RKKY interactions 
within the electronic Griffiths phase scenario. 

Our main results are as follows: (a) for disorder W 
weaker than a critical value W c we find Fermi liquid be- 
havior, but for W > W c the Kondo effect is suppressed 



on a finite fraction of spins, resulting in a spin- liquid 
phase displaying (universal) marginal Fermi liquid be- 
havior; (b) the spins that remain screened are still char- 
acterized by a power law distribution of (renormalized) 
energy scales P (T*) — (7 1 *)" ) but the exponent a* 
acquires a universal value a* « 1/2 throughout the spin- 
liquid phase; (c) the spin liquid phase is unstable to spin 
glass ordering at the lowest temperatures, but we find 
robust marginal Fermi liquid behavior in a broad tem- 
perature window above the freezing temperature. 

We consider the disordered Kondo lattice model as 
given by the Hamiltonian 

H = - l E (4 c ia + H - C + E V i 4 C * 

(ij)<r io 

» («■) 

where the exchange couplings between localized spins, 
and site energies i>i are distributed according to Gaussian 
distributions, PjUij) ~ exp(— Jfj/2 J 2 ) and Pw( v i) ~ 
exp(— vf/2W 2 ) :6j. In this expression, Si and Sj = 
\ E a ^ 4 <T a/3 c j / a represent a localized spin and the con- 
duction electron spin density at site i, respectively. We 
concentrate on the paramagnetic phase. Applying the 
standard procedure to average over disorder in the Jy 
couplings |2J and taking the limit of infinite coordination 
z — ► oo |8|, the local effective action assumes the form 



A = E / dT i dT ' c Ur)[(d T - M + VjW r') 
a Jo Jo 

- t 2 G c (r - r')]c i(7 (r') + J K f drS^r) ■ s,(r) 

Jo 

~ Y J dT J dT '^ T r, )S,-(r) • S,(r'). (2) 



The local action of Eq. (2) describes the so-called Bose- 
Fermi Kondo model |3, |lfj, which describes a Kondo 
spin interacting with both a fcrmionic bath of conduction 
electrons and a bosonic bath of spin fluctuations. For a 
disordered Kondo lattice, we must consider an ensemble 
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of such impurity models supplemented by following self- 
consistency conditions. The bosonic spin bath x( T ) = 
(TVS, (r). 8,(0))^ = J dVjPwiv^Sjir) ■ S 3 (0)) Aj 
and the conduction electron bath G c (r) = G C j(r) — 

— {T t c- (J {t)c^ !J {Qi)) ^ are obtained by appropriate disor- 
der averaging [fj , and for simplicity we use a simple semi- 
circular model density of states for conduction electrons. 

Destruction of the Kondo effect. The presence of 
RKKY interactions introduces a qualitative modification 
in the dynamics of the Kondo spins, through the presence 
of a dissipative bosonic bath of spin fluctuations. This 
behavior depends crucially on the precise spectral form 
of the bosonic bath, allowing for the destruction of the 
Kondo effect in the presence of sub-Ohmic dissipation 
[pj, [l0( . For a spectrum of the form 

X(iuj n ) ~ x(0) - C\uj n \ l ~ e , 

Fermi liquid behavior is recovered for e = 0, but for 
e > (sub-Ohmic dissipation), and for sufficiently small 
bare Kondo temperature Tk, the spin decouples from 
the conduction electrons. Within an electronic Grif- 
fiths phase, however, the disordered Kondo lattice has 
a very broad distribution of local Kondo temperatures 
P(T K ) ~ (Tk)^ 1 . Therefore, for e > and arbitrar- 
ily weak coupling to the bosonic bath (i.e. weak RKKY 
interaction), a fraction of the spins will decouple. 

To obtain a sufficient condition for decoupling, we ex- 
amine the stability of the Fermi liquid solution, by con- 
sidering the limit of infinitesimal RKKY interactions. To 
leading order we replace x( T ) — ► Xo ( r ) = x( T i J = 0), 
and the calculation reduces to the "bare model" of 
Ref. [jj. The resulting bosonic bath, which is an aver- 
age over the site-dependent local dynamic spin suscepti- 
bility, Xo (iu n ) = J dT K P (T K ) x {Tk, iv n ) , has a Fermi 
liquid form in the presence of weak disorder. However, 
for stronger randomness, W > W* « y/t 2 p c JK/2 cor- 
responding to a < 2 (here, p c is the density of states 
for conduction electrons) , the power law distribu- 
tion of energy scales within a Griffiths phase produces 
sub-Ohmic dissipation, corresponding to e = 2 — a > 0. 
Note that the estimate based on the bare theory sets an 
upper bound for the true critical disorder strength, i.e. 
W c < W* = W nfl /V2 (here, W nf i « \A 2 p c Jk /2 is the 
threshold for NFL behavior in the bare model [6j, cor- 
responding to a — 1). We emphasize that within the 
electronic Griffiths phase, such decoupling emerges for 
W > W c even for arbitrarily small J, in contrast to the 
clean case 0] where much stronger RKKY interactions 
(J > J c w 10 Tk) are required to destroy the Kondo 
effect. 

The spin liquid phase. For finite J, the actual value of 
e has to be self-consistently determined, as follows. For 
W > W c , the spins break up into two groups: the decou- 
pled spins and those that remain Kondo screened. Since 
the self-consistent bosonic bath function x(ia> n ) is an al- 



gebraic average over all spins, it is an additive function 
of the contributions from each fluid 

X {ioj-n) = n Xdc (iv n ) + (1 - n)x s (iw n ) ■ (3) 

Here, n is the fraction of spins in the decoupled phase. 
As we shall see, the functions Xdc (iu n ) and Xs («w„) 
both have a singular, non-Fermi liquid form character- 
ized by exponents £d c and e s , respectively. Deferring for 
a moment the study of the critical region (infinitesimally 
small n), we first examine the solution deep within the 
spin liquid phase. The first step in the self-consistent 
procedure is computing Sdc and e s for a given value of 
the bath exponent e. The spin autocorrelation func- 
tion in the decoupled phase assumes the form [j| [^j 
Xddr) = (TVS (r) • S (0)) - 1/t £ , a result valid to all or- 
ders in e [f(|. Since £d c is defined by Xdci?) ~ l/T 2 ~ Edc , 
we find 

£ dc (e) =2-e. (4) 

The non- analytic part of Xs {iw n ) comes from the spins 
with the smallest (renormalized) Kondo temperatures T* 
( "barely screened spins" ) 

Xta(»w»)= / dT*P(T*) Xbs {T*,icJ n ). (5) 
Jo 

Here P(T*) is the distribution of renormalized Kondo 
temperatures (local Fermi liquid coherence scales), and 
Xbs(T*,iuj n ) is the local dynamic susceptibility for a 
given T*. Properties of the Bose- Fermi Kondo model 
in the critical region of the decoupling transition have 
been extensively studied within renormalization group 
(RG) 0, and large- N approaches 01 1 and we use 
these results to calculate Xbs- In particular, T* ~ 
(5J K ) V - (STkY, which gives dT*/dT K - (T*) 1_1/l/ . 
Therefore P (T*) = P [T K (T*)] dT K /dT* - dT K /dT* ~ 
(T*) 17 ^ 1 . From scaling arguments HE3, Xbs (T*,uj) = 
(T*) 11 " 1 (j){u/T*), where 1] is the anomalous dimension, 
which is known to be exactly e ^(J ■ Performing the inte- 
gration in Eq. we find at low frequencies Xbs («w„) = 
Xbs(0) — C" \oj n \ v+ "~ , or, equivalently, at large times 
Xf«( T ) ~ l/r v+ i . By definition, Xs( T ) ~ l/ r2_£s ! which 
gives 

e s (e) = 2-r]-l/u. (6) 

Since rj — e and v > (as the relevant eigenvalue at 
the unstable fixed point), Eqs. (@J and © imply that 
£dc > £«< Therefore, £ = max{ed c , e s } = £dc, and from 
Eq. Q we find that the self-consistent bath is character- 
ized by the exponent e = 1, as in the spin liquid model of 
Sachdev and Ye , producing a logarithmic divergence 
of the average local dynamic susceptibility. Note that, in 
contrast to the bare (J = 0) model of the electronic Grif- 
fiths phase, the renormalized distribution P (T*) of local 
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energy scales now assumes a universal form character- 
ized by an exponent a*(e) = l/v(e) w e/2 = 1/2 within 
the spin liquid phase. More work is needed to determine 
the behavior of the uniform susceptibility, as well as the 
behavior of the specific heat. 

Transport in the spin liquid phase. Although the renor- 
malized Kondo coupling scales to zero for the decoupled 
spins, the precise form of the RG flows (scaling dimen- 
sion of "irrelevant operators") near the spin- liquid fixed 
point still determines the finite frequency (or finite tem- 
perature) corrections. To leading order, the contribution 
from decoupled spins scales as p(co) ~ [J k (lj)) 2 , while 
J K {to) ~ . To compute the appropriate exponent at 
the spin-liquid fixed point we have used the e-expansion 
approach of Ref . fip]| , and we find v = 2/e + 0(e 3 ). 
From our self-consistent solution for the spin-liquid phase 
(e = 1), we obtain \jv w 1/2, producing again a marginal 
Fermi liquid correction to the resistivity 5pd c (u) ~ oj, or 
at u> — and finite temperature 

Spdc(T) ~ T. 

Numerical results. As an illustration of our analytical 
predictions, and to obtain quantitative results, we pro- 
ceed to the numerical solution of our equations in the 
large- TV limit [ill fl3| . Introducing site-dependent slave 
boson parameters Tj and Sfj, and minimizing the local 
free e nerg y, we come to the following saddle-point equa- 
tions [nf 
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Figure 1: Local dynamic magnetic susceptibility, (a) At low 
frequencies, the decoupled spins (dashed line) provide the 
leading logarithmic dependence of the total averaged sus- 
ceptibility (full line). For the values of the parameters used 
(Jk = 0.8, J = 0.05, fi = —0.1 in units of the half bandwidth, 
corresponding to Tk {vj =0) =0.1 and W c « 0.1), there are 
n = 8% of decoupled spins at W — 0.4. The bare model 
(J — 0) leads to a stronger non- universal power law singu- 
larity (dash-dotted line), (b) x(* w ) f° r the disorder strength 
ranging from to 0.4. 



Critical behavior. Near the critical point the ar- 
guments which followed Eq. have to be modified 
since the relative importance of the various contribu- 
tions to the average local susceptibility changes. First, 
we concentrate on the contribution from the barely 
screened spins given by Eq. (J5J. As before, P (T*) = 
P [T K (T*)] dT K /dT*, but close to the transition P (T K ) 
is small and cannot be replaced by a constant prefactor 
of order 1. Since P (T K ) w P (T Kc ) ~ (T^)^ 1 , we find 



The local /-pseudo-fermion Green's function G fj {r) = « = Jo dT K PjT K ) ~jT K? ) a , where T Kc is the bare 
-{T T fj a {T)f] a (0)), is given by G7 1 (iuj n ) = iu n - 



(W n ) 



|A/j (iuj n ). The self-energy is equal 
to Ej-(r) = J 2 x(T)G fj (T), and A~f(iuj n ) = iu> n + 



/j, — Vj — t 2 G c {iu> n ). Self-consistency requires x( T ) — 
and G c (iLo n ) = G c j(iu) n ), where 



-Gn(r)G n (-r), 

G cj (iuj n ) = A f j (iuj n ) - r 2 / [iuj n - e fj - (iuj n )} . 

These equations were solved on the imaginary axis at 
T = using fast Fourier transform methods. The total 
average local dynamic susceptibility \ together with the 
contributions coming from Kondo screened \s and decou- 
pled spins \dc is shown in Fig. [TJa). At low frequencies, 
the contribution from Kondo screened spins saturates to 
a constant, while the decoupled spins produce a logarith- 
mic divergence. A comparison with the bare model illus- 
trates how the strong power law divergence of x found 
for J = is suppressed by the dynamical RKKY interac- 
tions. Fig. n^b) shows how x evolves with the change of 
disorder. Note that marginal Fermi liquid behavior per- 
sists up to a crossover scale lo s i ~ Q.1Tk{Vj = 0) which 
has very weak dependence on the disorder strength. 



Kondo temperature at the site energy v c at which the 
spins start to decouple. Therefore P {Tk) ~ n( a ~ x ^ a . 
From the bare model, we know that (for small J) a » 2 
near the critical point. Now we are in a position to write 
down the general form of the total bosonic bath at low 
frequencies 



X(iu n ) =Xo-Ci \uj„\ - C 2 n 1/2 \co n \ 



C 3 n In |cj n | . 

(9) 



The first two terms come from the well screened spins 
and have the Fermi liquid form. The third term is due 
to the "barely screened" spins and the last term is the 
contribution from the decoupled spins. The crucial point 
is that the non-analytic term from the barely screened 
spins, being proportional to ^/n, is much larger than the 
logarithmic term due to the n decoupled spins, except 
at exponentially small frequencies. Therefore, we can 
neglect the last term in Eq. 0. Below the crossover 
frequency u>* ~ n ( 1 / 2 )/[ 2 - i ;-( 1 / l/ )] j the non-analytic term 
in Eq. @ is dominant and close to the decoupling point, 
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Figure 2: Phase diagram obtained for the same values of pa- 
rameters as in Fig. The inset shows the fraction of decou- 
pled spins as a function of disorder. Note that the decoupling 
sets in already in the presence of moderate disorder. 



i.e. in the limit e/j — > and rj — > 0, Gf assumes the form 
Gf(iu>) « -iCsgncj/|cj| 1_£/2 for \u>\ < u>* , and Gf(iu>) ps 
— i/ui for \u>\ > u>*. The parameter e has to be self- 
consistently determined from the equation 1 — e = r\ + 
~ — 1. Within the large- N theory, v — 2/e which gives e = 
4/5. Inserting this expression into Eq. |JSJ, we find the 
critical site energy for decoupling v c ~ y| lnn|. Since the 
number of decoupled spins is equal to n = J v dvPw(v), 
we have a closed set of equations for n(W), from which we 
find that the number of decoupled spins is exponentially 
small in the vicinity of the critical point 



-A/{W-W C ) 



(10) 



where A is a positive constant. The numerical results for 
n(W) are shown in the inset of Fig. Interestingly, no 
precursors (vanishing coherence scale of the Fermi liquid) 
arise as the critical point is approached from the FL side, 
in contrast to what we have found by solving the same 
equations in absence of disorder 0]. This indicates a 
novel type of quantum critical behavior that has a char- 
acter of an essential singularity, a feature that appears 
specific to quantum Griffiths phases. 

Spin glass instability and phase diagram. In this pa- 
per we have concentrated on the paramagnetic solution 
of our model. However, the decoupled spins can be ex- 
pected to form a spin glass (SG) at low temperatures in 
the presence of random inter-site interactions [3. For a 
rigorous treatment of the spin glass phase, one needs to 
go beyond the N — oo limit, but a rough estimate of 
the temperature for SG ordering [JjJ may be obtained 
by using the large- N approach as an approximate theory 
for the considered N = 2 case. The spin glass instability 
criterion 0, as appropriately generalized to the case of 
additional site randomness then reads 



Fig. |3 represents a generic phase diagram of our model. 
For weak disorder the system is in the Fermi liquid phase, 
while for W > W c the marginal Fermi liquid phase 
emerges. The crossover temperature (dashed line) de- 
limiting this regime can be estimated from the frequency 
up to which the logarithmic behavior in x(icj) is observed 
UJ. The spin glass phase, obtained from Eq. \1 lfl . ap- 
pears only at the lowest temperatures, well below the 
marginal Fermi liquid boundary |l4j |. Interestingly, re- 
cent experiments have indeed found evidence of dynami- 
cal spin freezing in the milliKelvin temperature range for 
some Kondo alloys (lif . 

To summarize, we have introduced and solved a disor- 
dered Kondo lattice model with random inter-site RKKY 
interactions. Our solution, valid within extended dynam- 
ical mean-field theory, illustrates how non-Ohmic dissi- 
pation arising from inter-site RKKY interaction restores 
universality for non-Fermi liquid behavior of electronic 
Griffiths phases. Although considerably different in de- 
tail, this dissipative mechanism is reminiscent of the pro- 
cesses leading to dynamical freezing of droplets within 
magnetic Griffiths phases 0, suggesting a generic role 
of RKKY interactions in disordered heavy fermion sys- 
tems. 
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